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Abstract 

We study asymptotic synchronization at the level of global attractors in a class of coupled 
second order in time models which arises in dissipative wave and elastic structure dynamics. 

Under some conditions we prove that this synchronization arises in the infinite coupling inten¬ 
sity limit and show that for identical subsystems this phenomenon appears for finite intensities. 

Our argument involves a method based on “compensated” compactness and quasi-stability 
estimates. As an application we consider the nonlinear Kirchhoff, Karman and Berger plate 
models with different types of boundary conditions. Our results can be also applied to the non¬ 
linear wave equations in an arbitrary dimension. We consider synchronization in sine-Gordon 
type models which describes distributed Josephson junctions. 
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1 Introduction 

Our goal in this paper is to study long-time dynamics of a class of coupled systems consisting 
of two second order in time evolution equations. These systems are abstract models for studies 
of elastic and wave dissipative dynamics in different situations. Under some set of hypotheses 
concerning the model we first prove the existence of finite dimensional global attractors and study 
their dependence on interaction operators. Then we apply these results to analyze synchronization 
phenomena. In this paper we understand these phenomena at the level of global attractors. This 
means that in synchronized regime the attractor of coupled system becomes “diagonal” in some 
sense. 

Our main results are presented in Theorems 3.12 and 3.14. In particular, Theorem 3.12 
proves asymptotic synchronization in the limit of large coupling and Theorem 3.14 dealing with 
interaction of identical systems shows that synchronization is possible for finite values of cou¬ 
pling intensity parameter. As a preliminary step we obtain a result on uniform dissipativity with 
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respect to coupling intensity parameters (see Theorem 3.5). We also discuss possibility of syn¬ 
chronization in infinite-dimensional systems by means of finite-dimensional interaction operators. 
As application of these results we consider a range of nonlinear elastic plate models and also wave 
dynamics of different types. 

We note that recently the subject of synchronization of coupled (identical or not) systems has 
received considerable attention. There are now quite a few monographs [4, 28, 30, 34, 38, 40] in 
this field, which contain extensive lists of references. In the case of infinite dimensional systems 
the synchronization problem has been studied in [5, 7, 8, 24, 36] for coupled parabolic systems. 
Synchronization in Berger plates (they are a particular case of our abstract models) was considered 
in [31, 32, 33]. Master-slave synchronization of coupled parabolic-hyperbolic PDE systems was 
considered in [10, 11]. The methods involved in these publications relies either on the parabolic 
regularization (see [5, 8, 24]) or on a special structure of the model (see [31, 32, 33]). Our approach 
is different. 

par As an important technical tool we involve the method developed in [19] (see also also [20, 
Chapter 8] and [12]) based on an observability-type estimate which allows us to establish uniform 
quasi-stability estimates in the case of critical nonlinearities (such as in the von Karman and 
Berger models). In the standard way (see, e.g., [12], [19] or [20, Chapter 8]) these quasi-stability 
estimates lead to appropriate uniform bounds for attractors which are important for asymptotic 
synchronization. 

The paper is organized as follows. The next Section 2 is devoted to preliminary considerations. 
We describe here our abstract model, formulate main hypotheses and provide well-posedness result 
for rather general situation. Section 3 contains our main results on attractors and synchronization. 
In Section 4 we discuss possible applications. The Appendix contains some general facts from 
the dissipative dynamical systems theory. 

2 Preliminaries 

In this section we describe the problem and state our basic notations and hypotheses. Then we 
provide a well-posedness theorem which is easily derived from the known results. 

2.1 Abstract model and main hypotheses 

In a Hilbert space H we deal with following coupled equations 

u u + viAu + D\iu t + D V2 v t + Kuu + K i2 v + Bi(u) = 0, 

v tt + v 2 Av + D 2 iu t + D 22 v t + K 2 \u + K 22 v + B 2 (v) = 0, 
equipped with initial data 

it(0) = n 0 , u t { 0) = m, v(0) = vo, v t (0) = v lt 

under the following set of hypotheses. 

Assumption 2.1 (i) A is a self-adjoint positive operator densely defined on a domain T>(A) 

in a separable Hilbert space H , v\, v 2 >0 are parameters. We assume that the resolvent 


(la) 

(lb) 

(lc) 
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of A is compact in H. This implies that there is orthonormal basis {e^} in H consisting of 
the eigenvectors of the operator A : 


Ae/c = Xk^k, 0 < Ai < A 2 < • • • , lim A k = 00 . 

k —^00 

We denote by || • || and (-, •) the norm and the scalar product in H. We also denote by H s 
(with s > 0) the domain V(A S ) equipped with the graph norm || • || s = ||A S • ||. H~ s denotes 
the completion of H with respect to the norm || • ||_ s = ||A _S • ||. The symbol (-, •) denotes 
not only the scalar product but also the duality between H s and H~ s . Below we also use 
the notation H s = H s x H s . 

(ii) The damping operators D t j : H 1//2 1 —> H~ 1 / 2 are linear mappings such that the matrix 
operator 

V = : H 1/2 x H 1/2 1 —y H~ 1/2 x H~ 1/2 

\L>2i JJ 21 J 

generates a symmetric nonnegative bilinear form on H 1/72 = H 1 / 2 x H 1 / 2 . 

(iii) The interaction operators Kij : H 1 ^ 2 1 —> H are linear and 

lc= (K } i M . x H 2 ! 2 1 — y H x H 

\JX 21 A 21 / 

generates a symmetric nonnegative bilinear form on H 1 / 2 . 


(iv) 


The nonlinear operators Bi : H 1//2 — > H are locally Lipschitz, i.e., for every q > 0 there 
exists a constant L(o) such that 


II Bi(u) - Bi(v) || < L{q)\\u — u||i/ 2 , i = 1,2, 

for all u, v £ H 1 / 2 such that |MIi/ 2 ) IMI 1/2 — Q- addition we assume that Bi(u) are 

potential operators. This means that B t (u) = n((u), where Ilj : H 1 / 2 1 —> M is a Frechet 
differentiable functional on H 1 / 2 , i.e., 


Ihl 


lim 

1/2 


->o m 


II 1/2 


Bi(u + v) - Ilj(u) - (n'(u), v) 


= 0 . 


We also assume that II,(y) = noi('u)+IIij('u), where Iloj(u) is a nonnegative locally bounded 
functional on H l t 2 and 


\/ij > 0 3Crj : |n u (u)| < r/[||A 1 / 2 u|| 2 + U 0i (u)] + u^H 1 / 2 . 

The problem in (1) can be written as 

U tt + AU + VU t + K.U + B(U) = 0, U(Q) = U 0 , U t (0) = U u (2) 


where the operators T> and /C are defined above and 


U = 



A 





As it was already mentioned in Introduction our main motivation for (1) (and (2)) and also for 
the hypotheses in Assumption 2.1 is related with applications to coupled plate and wave systems 
(see Section 4). 
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2.2 Well-posedness 

Now we consider the existence and uniqueness theorem for (2). We start with adaptation of the 
standard definition (see, e.g, [19] and the references therein) to our model. 

Definition 2.2 A function [7(f) G C([0, T]; 77 1 / 2 )nC' 1 ([0, T]; 77) possessing the properties U(0) = 
Uq and Ut( 0) = U\ is said to be 

(S) strong solution to problem (2) on the interval [0, T], iff 

• u £ Wi(a, b ; H 1//2 ) and ut G W\ (a, 6; 77) for any 0 < a < b < T, where 

WiV, &; #) = {/ G C(a, 6; H) : f'€ Li(a, b;H )} , 

• AU(t ) + DUt(t ) G H for almost all t G [0,T]; 

• equation in (2) is satisfied in 77 for almost all t G [0,T]; 

(G) generalized solution to problem (2) on the interval [0, T], iff there exists a sequence of strong 
solutions {U n (t)} with initial data ( Uo n ,Ui n ) instead of (t/o, U\ ) such that 

lirn max { \\d t U(i) - d t U n (t )|| + \\A 1/2 {U(t) - t/ n (f))||j = 0. 
n-ux>t e [o,T] f > 

Application of Theorem 1.5 from [19] gives the following well-posedness result. 

Theorem 2.3 Let T > 0 be arbitrary. Under Assumption 2.1 the following statements hold. 

• Strong solutions: For every (t/o; U\) G 77 1 / 2 x 77 1 / 2 , such that AUq + VU\ G 77 there 
exists unique strong solution to problem (2) on the interval [0, T] such that 

(■ Ut-,Utt ) G ^^(O, T ; 77 1 / 2 x 77), U t G G r ([0,T); 77 1 / 2 ), 

l/i(GC r ([0,T);7t) and .4t/(t) + VU t (t) G C r ([0, T); 77), 

where we denote by C r the space of right continuous functions. This solution satisfies the 
energy relation 

£(U(t),Ut{t))+ [ (T>U t {r ), U t (r))dr = £(Uq, U\), (3) 

Jo 

where the energy £ is defined by the relation 

£(U 0 ,U 1 ) = £i(u 0 ,ui) +£ 2 (v 0 , Vl ) + £ int (u 0 ,v 0 ), 
with Uq = (uo;v 0 ), U\ = (uy,v i), 

£i{u 0 ,ui) = Ei(u 0 ,ui) + Ui(u 0 ) = i (||ui|| 2 + z/j||A 1/2 u 0 || 2 ) +IIj(n 0 ). 

and ^ 

£int{F Oj^o) — -(/CtTo, C/ 0 ). 
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• Generalized solutions: For every (U 0 \Ui) € H 1 / 2 x H there exists unique generalized 
solution. This solution possesses the property V l / 2 ut € £ 2(0 ,T\H) and satisfies the energy 
inequality 

£{U(t),U t (f)) + f {VU t {T),U t {j))dT<S(JJ Q ,U{). (4) 

Jo 

If U 1 and U 2 are generalized solutions with different initial data and Z = U 1 — U 2 , then 
\\Z t (t)\\ 2 + \\A l / 2 Z(t)\\ 2 + \\lC l / 2 Z(t)\\ 2 < (||Z t (0)|| 2 + ||A 1/2 Z(0)|| 2 + \\K}l 2 Z(t)f) e aRt 

provided \\U)(0)\\ 2 + ||A 1 / 2 £P(0)|| 2 + ||/C 1 / 2 t/ i (0)|| 2 < R 2 . 

By Theorem 2.3 problem (2) generates a dynamical system (H, St) in the space 

U = h 1/2 x H = H 1/2 x H l > 2 x H x H 
with evolution operator defined by the relation 

St(U Q -,lh) = (U(tfU t (t)), 

where U(t) is a generalized solution to problem (2). 

3 Global Attractors 

In this section we prove the existence of a global attractor for the dynamical system fH, St) 
and study its properties. Keeping in mind further application we impose additional hypotheses 
concerning the damping operator T> and the source term B. 

Assumption 3.1 Let Assumption 2.1 be in force and 

(i) V is strictly positive, i.e. exists cq > 0 such that 

{VW,W)>c 0 \\W\\ 2 , WeH 1/2 ; 

(ii) either Bt are subcritical, i.e., for every g > 0 there exists a constant L(g) such that 

\\Bi(u) - Bi(v)\\ < L(g)\\u-v\\ 1/2 _ s , i = 1,2, 6 > 0, 

for all u, v £ H 1 / 2 such that || , u|| 1 / 2 ! IMI 1/2 < f? 5 or else the potential energies II, are 
continuous on H 1 / 2- ' 5 for some <5 > 0 and the mapping u 1 —> A~ l Bi{u) is continuous from 
H l / 2 ~ s into H~ l for some 6, l > 0, i = 1, 2. 

Proposition 3.2 Let Assumption 3.1 be in force. Then the system {Ft, St) generated by problem 
(2) is asymptotically smooth (see the definition in the Appendix). 

Proof. We apply Theorem 3.26 and Proposition 3.36 from [19]. Since 

|(W,W)| <[{VV,V)] 1/2 [(VW,W)} 1/2 < C[(VV,V)] 1/2 \\A 1/2 W\\ 

<C £ {VV,V) + e\\A 1/2 W\\ 2 , 

relation (3.60) in [19, p.54] obviously holds in a simplified form. □ 
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One can see from the energy inequality in (4) that the system (LI, St) is gradient with the 
full energy £(Uq\U\) as a strict Lyapunov function (see the Appendix for the corresponding 
definitions). Therefore by Corollary 2.29 [19] (see Theorem A.2 and Remark A.3 in the Appendix) 
to guarantee the existence of a global attractor we need boundedness of equilibria. This leads to 
the following assertion. 

Theorem 3.3 Let Assumptions 3.1 be in force. Assume in addition that there exist u < Ui and 
C > 0 such that 

z/||A 1 / 2 'u|| 2 + u) + C > 0 , u € H 1 ^ 2 , i = 1,2. (5) 

Then the system ( LL,St ) generated by problem (2) possesses a compact global attractor. 

Proof. Stationary solutions U = (it; v) € H l t 2 solve the problem 

v\Au + I\\ 2 U + K\ 2 V + B\(u) = 0, 
u 2 Av + K 2 \u + K 22 v + B 2 (v) = 0. 

Using the multipliers u for the first equation and v for the second and also positivity of the 
operator 1C, we obtain that 

z^i||yl 1//2 tt|| 2 + ^ 2 ||^ 1/2 t’|| 2 + ( Bi{u),u) + (. Bi(v),v ) < 0 

By (5) this yields ||M 1 / 2 u|| 2 + ||M 1 / 2 u|| 2 < C (with C independent of V and 1C). Thus stationary 
solutions are bounded. □ 

3.1 Uniform dissipativity 

For synchronization phenomena it is important to have bounds for the attractor independent of 
interaction operators T> and 1C. In spite of the set of stationary solutions is uniformly bounded 
with respect to V and 1C Theorem 3.3 does not provide appropriate bounds for the attractor. 
Below we use an approach based Lyapunov type functions which allows us to prove uniform 
dissipativity of the system (LI, St). To simplify argument it is convenient to introduce intensity 
factors a and x for interactions in velocities and displacements. Moreover, we assume a particular 
structure of V related with interaction operator 1C. So instead of (2) we consider 

U u + AU + {V 0 + alQUt + xlCU + B(U) = 0, U(0) = U 0 , U t (0) = U u (6) 

where a and x are nonnegative parameters. 

In addition to Assumption 2.1 we impose the following hypotheses. 

Assumption 3.4 (i) The operator Vq is bounded from H 1 / 2 into H and there exist Co > 0 

and a > 0 such that 

({V 0 + aJC)W,W) > c 0 ||IU|| 2 , WeH 1 ' 2 - 

(ii) there exists r] < 1 and 6, C > 0 such that 1 

5no(U) - (B(U),U) < ri\\A 1/2 U\\ 2 + C , U^H 1 ' 2 , (7) 

where flo(?7) = Iloi(it) + IIo 2 (u) with JJ = {u; v). 

1 The relation in (7) is the standard requirement in many second order in time models, see, e.g., [19]. 
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Theorem 3.5 Let Assumptions 2.1 and 3-4 be in force. Then for every a > a and x > 0 
the system (Ti,St) generated by problem ( 6 ) is dissipative 2 with an absorbing ball of the size 
independent of (a; x) £ A = {a > a, x > 0 }. More precisely, there exists R independent of 
(a; x) G A such that the set 

<3 = {(U 0 ]U l )e'H : ||C/i || 2 + m 1 / 2 C/o || 2 + x||/C 1 / 2 C / 0 || 2 < i? 2 } ( 8 ) 

is absorbing. 

We note that the estimate in ( 8 ) improves the corresponding finite-dimensional statement in [ 1 ] 
which requires uniform boundedness of the ratio a/x. As it is shown in Proposition 3.8 we need 
the latter property for uniform quasi-stability only. 

Proof. We use a slight modification of the standard method (see, e.g., [3, 9, 39]) based on 
Lyapunov type functions. Let U(t) = (u(t)",v(t)) be a strong solution, 

E(t) = E(U ; U t ) = ^ (\\U t (t)\\ 2 + \\A l ' 2 U(t)\\ 2 ) + U(U(t)) 

with n([/) = Ili(tt) -TII 2 (i’), and <h(f) = r](U, Ut) + /i(1CU, U), where 77 is a positive constant which 
will be chosen later and 2/i = x + pa. We consider the functional V = E + <h. One can see that 
there exist 0 < 770 < 1 and /% > 0 independent of (a; x) such that 

A) [E*(t) + x||/C 1 // 2 C/(f)|| 2 ] -Pi<V< /3 2 \E*(t) + p\\JC 1/2 U(t)\\ 2 ] +/3 3 (9) 

for all i] € ( 0 , 770 ], where 

E*(t) =E,(U;U t ) = ^ (\\Ut(t)\\ 2 + ||^4 1 / 2 C/(t)|| 2 ) + n 0 (J7). 

Now on strong solutions using the energy relation in (3) we calculate the derivative 

^ = -((Vo + aJC)U t , Ut) + plWUtW 2 - (V 0 U t , U) - (AU, U) - x(KU , U) - (B(U), U)]. 

Since Vq is bounded from H 1 / 2 into H, we obtain that 

\(V 0 Ut, U )I < e\\A 1/2 Uf + Ce~ l \\Ut\\ 2 , Me > 0. 

Thus by Assumption 3.4(ii) there exist bt> 0 independent of (a, x) such that 

^ < -[((V 0 + alC)U t ,Ut) - hpWUtf] — 6277 [E* (t) + x(KU , U)\ +7763, 

This implies that there exists 0 < 77* < 770 independent of (a, x) G A such that 

^ + hr,[E.(t) + x(JCU, U)\ < 7763 ( 10 ) 

2 See the Appendix for the notions related with this theorem. 
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for all (a, x) G A and // G (0, r/*] , where I/g = V + fa > 0. 

Now we split the parametric region A into several subdomains. 

We start with the following case. Let x* > 0 and a, > a be fixed. We take 0 < x < x* such 
that a* > x/?y* Then we take q = xa _1 . In this case q < and 

1 , , 1 , 

/x = — (x + jya) = -(x+x) < x. 

Thus for all x > x* and a > a* we have that 

^ + 6 2 ?? [£*(f)+MKI/,[/)] <qb 3 . 

In particular, (9) yields 

< hfa with b 4 = b 3 + ^ ~ , 

at p 2 

where 6 ,; and fa do not depend on a and x. This implies that 

V^(t) < Vp(0)e~ , » r ’K lt + 

02 

Since the value bfaof 1 fa is independent of x > 0 and a > a, we can conclude from the previous 
argument that the set 21 is absorbing for all x > 0 and a > a with R independent of x and a. 
In the case when x = 0 from (10) we obtain 

+ b 2 rjEfat) < rjb 3 . 

Now we take a > and rj = a -1 . In this case 77 < ?y* and 77 = 1/2. Therefore the conclusion 
follows from the estimate ||/C 1 // 2 I/|| < c||^I 1//2 [/|| by the same argument as above. 

In the case x = 0 and a < a < rfa 1 the conclusion is obvious. 

So the remaining case is x > 0 and a = a. Now we can take ij = min{ 7 /*, xd" 1 } when a > 0. 
It is clear that fj, < x for this case. Thus we can argue as above. In the case a = 0 the relation 

H < x/2 holds automatically. This completes the proof Theorem 3.5. □ 

This theorem and Proposition 3.2 immediately imply the following result on the existence of 
a global attractor. 

Theorem 3.6 Let Assumptions 2.1, 3.1(H) and 3-4 be in force. Then for every (a; x) € A the 

system (' H , St) generated by problem ( 6 ) possesses a compact global attractor 21 . For every full 

trajectory Y = {(£/(£); Ufatj) : t G M} from the attractor 

sup {||t/s(t)|| 2 + ||.4 1 / 2 t/(f)|| 2 + x||/C 1/,2 t/(t)|| 2 j + [°° (V 0 + aK,)Ufar), Ufar)) dr < R 2 (11) 

te K J J -00 

for some R independent of (a; x) G A. 

Proof. We first apply the standard result on the existence of a global attractor, see Theorem 
A.l in the Appendix. This attractor belongs to the set 23 defined in ( 8 ). This implies an uniform 
bound for the supremum in (11). Using the energy relation in (4) we obtain the corresponding 
bound for the dissipation integral in ( 11 ). □ 
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3.2 Quasi-stability 

The uniform bounds for the attractor given by Theorem 3.6 are not sufficient to perform the large 
coupling limit a —>• +oo and/or x —>• +oo in the phase state of the system. One of achievements 
of this section is stronger uniform estimates for the attractor size. For this we apply the quasi¬ 
stability method in the form suggested in [19] (see also [13, 15, 16, 18, 20] and the survey in the 
recent monograph [12]). This method makes it also possible to proof finite-dimensionality of the 
attractor and obtain its smoothness properties. 

To apply the quasi-stability method we need additional hypotheses concerning the nonlinear 
forces Bi{u). 

Assumption 3.7 Assume that 

• Bi(u ) = n'(it) with the functional H : H 1 / 2 i —> R which is a Frechet C7 3 -mapping. 

• The second II) (u) and the third II) (u) Frechet derivatives of IT (u) satisfy the conditions 

(nfV);u,u)| < C p \\A°v\\ 2 , v £ FT 1/2 , (12) 

for some a < 1/2, and 

(nf^u);^,^,^) < C , p||A 1/2 ui||||A 1/2 U2||||u 3 ||, u; <E tf 1/2 , (13) 

for all u £ H 1//2 such that ||A 1 / 2 u|| < p, where p > 0 is arbitrary and C p is a positive 
constant. Here above (n^(u);ui,... ,Vk) denotes the value of the derivative n (u) on 
elements vi,... ,Vk- 

This assumption concerning nonlinear feedback forces B{ (u) appeared earlier in the case of systems 
with nonlinear damping (see [19, p. 98] and also [14]) to cover the case of critical nonlinearities. 
We note that Assumption 3.7 holds in both cases of the von Karman and Berger models (see [19] 
p. 156 and p. 160 respectively). Moreover, as it is shown in [19, p. 137], this Assumption 3.7 is 
also true in the case of the coupled 3D wave equation in a bounded domain 9 C R 3 of the form 

u tt + °iu t - Au + ki(u - v) + (pi(u) = u| an = 0, (14a) 

Vtt + (J 2 Vt - Av + k 2 (u - v) + tp 2 {u) = f 2 (x), w| an = 0, (14b) 

provided ipi £ C 2 (R) possesses the property |<^/(s)| < C(1 + |s|) for all s £ R, the parameters 

ai and ki are nonnegative. Thus our abstract model covers the case of 3D wave dynamics with 
a critical force term. We refer to Section 4.2 for a further discussion concerning nonlinear wave 
equations. 

Recall that the Frechet derivatives (u) of the functional n are symmetric /c-linear continu¬ 
ous forms on H 1//2 (see, e.g., [6]). Moreover, if n £ C* 3 , then ( B(u),v ) = (H' (u)\v) is C 2 -functional 
for every fixed v £ H 1//2 and the following Taylor’s expansion holds 

(B(u + w) — F(u), v) = (n^(u); w, v) + f (1 — (u + Xw); w, w, v)d\ (15) 

Jo 
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for any u, v € H 1 / 2 [6]. If Assume that u(t) and z(t) belong to the class C 1 (a, b; V(A 1 ^ 2 )) for 
some interval [a, 5] C M. Then, by the differentiation rule for composition of mappings [6] and 
using the symmetry of the form II^ (u) , we have that 

^<n (2 \u);z,z) = {U (3 \u)-,u t ,z,z) +2(Il {2 \u);z,zt). 

Therefore from (15) we obtain the following representation which is important in our further 
considerations: 

(Bi(u(t) + z(t)) - Bi(u(t)),zt(t)) = ^Qi(i) + Ri{t)r te[a,b\ CM, (16) 

with 

Qi(t) = -^(n^ («(*)); *(*)»*(*)> ( 17 ) 

and 

Ri(t) = --(Jlf\u)\ut,z,z) + / (1 - X)(uf\u +\z)-z,z,zt)d\. (18) 

1 Jo 

As we will see below the representation in (16) and the hypotheses listed in Assumption 3.7 can 
be avoided if the nonlinear forces Bi(u) are subcritical, i.e., 

3oo< 1/2: \\Bi{ui) - Bi{u 2 )\ < L(g)\\A ao (ui - u 2 )\\, V||^l 1/2 Uj|| < g. (19) 

The representation in (16) leads to the following assertion which, in fact, is proved in [19] (see 
(4.38), p. 99), but without control of the parameters a and k. 

Proposition 3.8 Let Assumptions 2.1 and 3-4 be in force. In addition we suppose that either 
Assumption 3.7 or relation (19) holds. Let A4 C B be a bounded forward invariant set with 
respect to St and U l = StU (, i = 1,2 be two solutions to (6) with (different) initial data U( G A4. 
Let Z = U 1 — U 2 . Then there exist C, 7 > 0 such that 

E z (t ) < CE Z ( 0)e" 7t + Cm&x\\A a Z(T)\\ 2 , Vt > 0, (20) 

[04] 

where 0 < cr < 1/2 and 

E z (t) = \ [\\Z t (t)\\ 2 + \\A l ' 2 Z{t)\\ 2 + x\\K}/ 2 Z\\ 2 ) . 

If M is a uniformly bounded in TL with respect to (a;x) € A and 

(a; x) € A^ = {(a; x) € A : a < (3(1 + x)} 

for some (3 > 0, then the constants C, 7 independent of (a; x) but can depend on (3. 

Proof. We use the same line of argument as in [19] and start with the following relation (which 
follows from Lemma 3.23 in [19]): 

E z (t)dt <c^yj ((A) + alC)Z t ,Z t )dt 

+ / \((D 0 + aIC)Z t ,Z)\dt + ^ T (U\U 2 ) 

Jo 


TE Z (T) + [ T 
Jo 
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for every T > T 0 > 1, where c > 0 does not depend on a, x, T and 



r T 


r T 


r T r T 

*t(U\U 2 ) = 

/ (G(r),Z t (r))dT 
Jo 

+ 

/ o G(t),Z(t))dt 
Jo 

+ 

dt Z t {r))dr 

Jo Jt 


with 

G{t) = B{U\t))-B(U 2 (t)). 

Since every point (a; x) G A belongs to for some (5 > 0. it is sufficient to consider the case 
when (a; x) € for some /?. In this case we have that 

|((A) + aiqz t , Z) I <((D 0 + a/C)Z t , Z t ) l / 2 ((D 0 + a/C)Z, Z) 1 / 2 

+ alC)Zt , Zt) + sEz(t ) 


for every e > 0. Thus choosing e in an appropriate way we obtain that 

TE Z (T)+ [ E z (t)dt < op [ {{D 0 + aK)Z t ,Z t )dt + c 0 ^ T {U\U 2 ). 

Jo Jo 

Under Assumption 3.7 we have from Proposition 4.13 in [19] that for any e > 0 and T > 0 there 
exist a(e,T) = a^(e,T) and b(e,T) = bjvi(£,T) such that 


sup 

te[o,T] 


(G(r),Z t (r))c?r 


rT r 


<e 


|Z,(r)l| 2 + l|V/ 2 Z(r) 


dr 


( 21 ) 


+ T ) 



d(r)||Al 1 / 2 Z(r)|| 2 dr + b(e,T) sup ||.A ct Z(t) || 2 

re[0,T] 


for all s > 0, where a < 1/2 and 


d{t)^d{t-U 1 ,U 2 ) = \\U}{t)\\ 2 + \\U 2 {t)\\ 2 . 


Obviously the same relation (21) (even with a(e, T ) = 0) remains true in the subcritical case (19). 
Thus 

y T (U\U 2 ) <e [ T E z (T)dT + a(£,T) l T d(r)P 1 / 2 Z(r)|| 2 dT + 6(e,T) sup M CT Z(r)|| 2 
Jo Jo re[0,T] 


for every e > 0. From the energy relation we also obtain that 

((A) + aK)Z u Z t )dt < E z { 0) - E Z (T) + V T (U\U 2 ) 
Thus after appropriate choice of £ and T we arrive at the relation 



E Z {T) < qE z (0) + a(T) 



d(r)m 1//2 Z(r)|| 2 dr + b(T) sup \\A a Z(T)\\ 2 , 

re[0,T] 


where q < 1 and all constants depend on fd. This inequality allows us to apply the same procedure 
as in [19, p.100] to obtain (20). □ 
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Now we are in position to obtain a result on the finiteness of fractal dimension of the attractors 
and also additional bounds for trajectories from these attractors. 

Theorem 3.9 Let Assumptions 2.1 and 3-4 be in force. In addition we suppose that either 
Assumption 3.7 or relation (19) holds. Then for any (a; x) G A the following assertions hold: 

1. The global attractor 2l" ,x of the system (fH, St) generated by (6) has a finite fractal dimension 
dim f 21. 

2. This attractor 2l a,x lies in H 1 ’xH l/2 and for every full trajectory Y = {(U(t);Ut(t)) : t G M} 
from the attractor in addition to the bound in (11) we have that 

sup{||C/ u (i)|| 2 + \\A l / 2 U t (t)f + x||/C 1 / 2 [/f(f)|| 2 j < Rti/3) (22) 

te R 1 J 

for some R\(/3) independent of (a; x) G A p, where f3 > 0 can be arbitrary. 

3. The attractors 2l" ,x are upper semicontinuous at every point (a*;x*) G A, i.e., 

lirn ["sup |dist w (y,2t"*’ x ‘) : y € VI 0 "’* 71 }} = 0 (23) 

n—foo 1 ^ J J 

for every sequence {(cd 1 ; x n )} C A such that (a n ; x n ) —> (a*; x*) G A as n —> oo. 

Proof. By Proposition 3.8 the system (H, St) is quasi-stable on every bounded forward invariant 
set. Thus we can apply Theorems 3.4.18 and 3.4.19 from [12] to prove the statements 1 and 2 
(see also Theorems A.4 and A.5 in the Appendix). 

To prove upper semicontinuity property we can use the methods developed in [25, 26]. In 
particular, we can apply a result due to [26] (see Theorem A.6 in the Appendix). Indeed, let 
( a n ; x n ) -» (a*; x*) G A as n —> oo. In this case it follows from (22) that 

\\AU(t)\\ 2 + \\A 1/2 U t (t)f < C 2 , Vt G K, 

for every full trajectory ( U(t ); Ut(t)) from the attractor 2l“ n,Xn , where C does not depend on n. 
This means that the attractor belongs to the set 

{(U 0 ,U i) : \\AUo\\ 2 + \\A l ' 2 U l \\ 2 <C 2 } 

which is compact in R. Thus we only need to show the property (ii) in Theorem A.6. 

Let {U(f\ U r () G 2t Q,l ’ Xn and (U 0 : U i) G 2t Q *’ x *. One can see that 

(Z(t)-z t (t)) = s? n ’*"(uz-Uf) - s“* ,x *([/o; C7i) = (U n (t) - u{t)- u?(t) - u t (t)) 

satisfies the equation 

Ztt + AZ + ( T)q + a*IC)Zt + xRQZ = F, 

where 

F = -K - a^ICU? - {x n - x*)lCU n - B(U n ) + B(U). 

On the attractors we obviously have that 

ll^ll < ci(|a n - a*| + \x n - x*|) + c 2 \\A l / 2 Z\\. 
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Therefore the standard energy type calculations gives the estimate 

|| S? n ’* n Y 0 n - S?*'"*Y 0 \\ n < C(\a n - a*| + |x n - x*| + ||T 0 ” - Y 0 \\ n )e at , t > 0, 

where Yq 1 = (Uq \ U f) and Y 0 = (C/ 0 ; U\). Thus we can apply Theorem A.6. □ 

Remark 3.10 The results stated in Theorems 2.3 and 3.9 deal with a general model of the form 
(2) or (6) and thus they can be also applied in the case of several interacting second order in time 
equations of the form 

N N 

u\ t + UiAu 1 + D 0i u l t + a ^ Kiju{ + x^ KijV? + Bi(u l ) = 0, i = 1,... ,N, (24) 

3 =1 3 =1 

under obvious changes in the set of hypotheses concerning the operators in (24). 

3.3 Asymptotic synchronization 

Now we apply the results above to synchronization. We switch on the interaction operators 1C of 
the standard (see, e.g., [8, 24, 31, 36] and the references therein) symmetric form. Moreover we 
suppose that the damping operator Vq has a diagonal structure. Namely we assume that 

K = (-1 4 1 ) = ( D o‘ d 2 ) < 25 > 

where K is a strictly positive operator in H with domain T>(K) D H 1 / 2 , the operators D{ : 
H 1 / 2 i —y H are nonnegative. Thus we consider the following problem 

utt + v\Au + D\Ut + aK(ut — vt) + xK(u — v) + B\{u) = 0, (26a) 

v tt + 1^2 Av + D 2 v t + aK(v t — u t ) + xK(y — u) + B 2 {y) = 0, (26b) 

with initial data 

it(0) = n 0 , u t { 0) = ui, n(0) = v 0 , v t (0) = v x . (26c) 

All theorems stated above can be applied to this situation. 

Our goal is to study asymptotic synchronization phenomena and we are interested in qualita¬ 
tive behavior of the system in the large coupling limit x —>• oo (and/or a —>■ oo). It is clear from 

the bound of the attractor given in (11) that it is reasonable to assume that u = v in the this 

limit. Therefore we need to consider a limiting problem of the form 

w tt + vAw + Dw t + B(w) = 0, u>(0) = w 0 , w t (0) = w±. (27) 

where 

v= Vl ^ ^ > D = ^i D i+ d -2), B(w) = ^{Bi(w) + B 2 (w)). 

Obviously the argument above can be applied to system (27) provided the damping operator D 
is not degenerate. In fact we can easily prove the following assertion. 
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Proposition 3.11 Let Assumptions 2.1(i,iv), and 3.4(H) be in force. Let D : H 1//2 H be a 
strictly positive operator. In addition assume that either Assumption 3.7 or relation (19) is valid. 
Then problem (27) generates a dynamical system in the space H 1 ^ 2 x H possessing a compact 
global attractor of finite fractal dimension. This attractor is a bounded set in H 1 x H 1 / 2 . 

Below we show that the attractors 2l" ,x for problems (26) in some sense converge to the attractor 
of the limiting system (27) when x —> +oo. 

Theorem 3.12 Let Assumptions 2.1, 3.1(H) and 3-4 be in force with 1C and Vq of the form 
given in (25). Then for every (a; x) € A the system (TL,St) generated by problem (26) possesses 
a compact global attractor 2l a,x . For every full trajectory Y = {(U(t);Uf(t)) : t. G R} with 
U(t) = (u(t);v(t)) from the attractor 

C + OO 

sup \\\U t (t)\\ 2 + \\A 1/2 U{t)\\ 2 + x\\K 1/2 (u(t)-v(t))\\ 2 \ + a / \\U t (T)\\ 2 dr < R 2 (28) 

t€M ^ ' ,/—oo 

for some R independent of (a; x) £ A. 

In addition assume that either Assumption 3.7 or relation (19) holds. Then for any (a; x) G A 
the following assertions hold: 

1. The global attractor 2l“ ,x of the system (' H,St ) generated by (26) has a finite fractal dimen¬ 
sion dimj2l. 

2. This attractor 2l a,x lies in H 1 xH 1 / 2 and for every full trajectory Y = {(U(t);Ut(t)) : t € K} 
from the attractor in addition to the bound in (28) we also have that 

sup { \\U t t{t)\\ 2 + ||^l i/2 t/ t (t)|| 2 + x\\K l > 2 (u t (t) - u f (t))|| 2 ) < R\{(3) (29) 

te R 1 J 

for some R\{j3) independent of (a; x) € Ag, where fi > 0 can be arbitrary. 

3. The attractors 2l" ,x are upper semicontinuous at every point (a*; xfi) € A, i.e., (23) is valid 

for every sequence {(a n ; x n )} C A such that ( a n ; x n ) (a*; xfi) € A as n —> oo. 

4. Let u i->- Bjf u) be weakly continuous from H 1 / 2 into some space H~ l , l > 0. Then in the 
limit x —>• oo we have that 

lim [sup {dist He (Y, 21) : Y G 2l“’ x jl = 0, (30) 

yc^r-oo L k _ 

where TL e = H x l 2 ~ £ x H l / 2 ~ £ and 21 = { (uo; uol ui', u\) '■ ( u o ; u i) € 2l}. Here 21 is the 
global attractor for the dynamical system generated by (27). Moreover, if instead of weak 
continuity of B, we assume that v\ = and K commutes 3 with A, then the convergence in 
(30) holds in the space H 1 / 2 x H x ^ 2 ~ £ C TL. 

3 We can take K = A a with some 0 < a < 1/2, for instance. 
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Proof. All results except the last one easily follows from Theorem 3.9. Thus we need to establish 
property (30) only. As in [25, 26] we apply contradiction argument. 

Assume that (30) is not true. Then there exist sequences {x n —>■ oo} and T 0 n € 2l" ,x " such 
that 

dist We (Y 0 r \2l) > <5 > 0, n = l,2,... 

Since T 0 n € 2P ,Xn , there exists a full trajectory Y n = {(C/ n (t); t/”(t)) : t G M} from the attractor 
2l a,Xn such that Y n ( 0) = lg n . It follows from (28) and (29) and also from Aubin-Dubinsky-Lions 
theorem (see [37], Corollary 4) that the sequence Y n is compact in C{a,b\'H e ) for every a < b 
and *-weakly compact in L 00 (M;77 1 / 2 x 77 1 / 2 ).Thus there exists 

U(t) = ( u(t ); v(t)) € C^K; U £ ) n L 00 (R-, H 1/2 x H l/2 ) with U t e L 00 (R- 1 H 1/2 x H 1/2 ) 

such that along a subsequence 

Va<6: sup {\\U?(t)-U t {t)\\\ /2 _ £ + \\U n (t)-U{t)\\\ /2 _\ -»■ 0, n ^ oo, (31) 

t£[a,b] 1 J 

and 

(U n (ty,U?(t)) -»■ *-weakly in L^H 1 ' 2 x H 1 / 2 ) , n —» oo. 

Since K is strictly positive, it follows from (28) and (29) that 

sup {IK(t) - v?(t)\\ 2 + I \u n {t) - v n (t) II 2 } < — 0, n oo. 

teR 

By interpolation 

sup - »,”(())« < Csup { (||A 1 / 2 <(«)II + l|. 4 1 / 2 <(«)ll) 1-2 ' K”(t) - 1>?(«)II 2 ' 

t€R t€R l v ' 

< Csup ||«”(t) — i’"(f)|| 2e —» 0, n -> oo. 
teR 

Similarly 

sup ||A 1,/2_£ (u n (t) — u n (t))|| —>• 0, n —> oo. 
te R 

Since u H>• Bi(u ) is weakly continuous for i = 1,2, these observations allow us to make a limit 
transition in the variational form of equations (26) and conclude that U(t) = (u(t)',u(t)), where 
u(t ) is a solution to (27). Moreover, (u(t);ut(t)) is a trajectory bounded in H 1 / 2 x Tf 1 / 2 . Thus it 
belongs to the attractor 21. It is also clear from (31) that 

Yo = (U n (0y, U?{ 0)) -»■ (U (0); ^(0)) € 2t in U £ 

which is impossible. 

In the case when v\ = v -2 and K commutes with A taking sum of equations (26a) and (26b) 
we can find that 

sup || A(u(t) + u(t))|| < C{Rp) (32) 

te R 
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for every trajectory (u(t);v(t)‘Ut(t);v t (t) from the attractor with (a; x) € A^. Taking the 
difference of (26a) and (26b) we obtain that 

sup ||Az(t) + 2xKz(t)\\ < C(R/3 ) with z(t) = u(t ) — v(t) (33) 

i€ R 

for fixed a. Since A and K commutes, 

|| Az + 2xKz\\ 2 = ||A;?|| 2 + Ax 2 \\Kz \\ 2 + 4 x(Az, Kz) 

= \\a z \\ 2 + ax 2 \\k z \\ 2 + ax\\a 1 / 2 k 1 i 2 z \\ 2 > pz || 2 . 

Thus (32) and (33) yield additional estimate on the attractor: 

sup{||Au(f)|| + ||An(t)||} < C 
te R 

with the constant C independent of x. This provides us with compactness of JJ n {t) in the space 
C([a, b\; H l ~ £ ) for every e > 0 and makes it possible to improve the statement in (30). □ 


Remark 3.13 • The result in (30) means that the attractor 2l“ ,x becomes “diagonal” in 

the limit of large intensity parameter x with fixed or even absent interaction in velocities. 
Thus the components of the system becomes synchronized in this limit at the level of 
global attractors. In particular, this implies that every solution U(t) = (u(t);v(t)) to (26) 
demonstrates the following synchronization phenomenon: 


Ve > 0 


lim sup 

t —^ ~|~oo 


- v t 




< £, V X > X* 


• The same conclusion as in (30) can be obtain in the limit (cr; x) —>• +oo inside of for some 
/3. However as one can see from (30) large a is not necessary for asymptotic synchronization. 
This observation improves the result established in [24] for finite-dimensional systems which 
requires for synchronization both parameters a and x to be large. 

• The possibility to obtain synchronization for fixed small x and large a is problematic. The 
point is that in the case x = 0 under appropriate requirement on nonlinear forces Bi there 
are possible two different stationary solutions which demonstrate absence of asymptotic 
synchronization. 

Now we consider the case of identical interacting subsystems, i.e., we assume that 

v\ = u 2 = v, Di=D 2 = D, Bi(w) = B 2 (w) = B{w). (34) 

In this case we observe asymptotic synchronization for finite values of x. 

Theorem 3.14 Let the hypotheses of Theorem 3.12 and also relations (34) be in force. Assume 
that a € [<5,a*] for some fixed cc*. Let 

s x = inf {i/(Aw,w) + x(Kw,w) : w € H 1//2 , ||u;|| = l} 
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There exists s* = s*(a, a*) and lo > 0 such that under the condition 4 s H > s* the property of 
asymptotic exponential synchronization holds, i.e., 

{ eUt II u t (t) ~ vt(t)\\ 2 + II A 1/2 u{t) - u(t))|| 2 |=0 (35) 


for every solution U(t) = (u(t)',v(t)) to (26). In this case 2t" ,x = 21 for all x such that s M > s*. 
Proof. In the case considered w = u — v satisfies the equation 

wtt + vAw + Dwt + 2 aKwt + 2 xKw + B(u) — B(v) = 0, w(0) = vjq, wt( 0) = vj\. 
where wq = uq — vq and w\ = u\ — v\. 

We consider the case of the critical nonlinearity, the subcritical case is much simpler. In the 
former case we use the representation (16) with z = w. Since B\ = £>2 = B, below we omit the 
subscript i. It follows from (12) and (13) that the variables Q and R defined in (17) and (18) 
admits the estimates 

1*2(01 < C R \\A°w(t) ||' 2 and |B(t)| < C R (||« t (f)j| + Mf)ll)IM 1/2 t»(t)|| 2 
under the condition 

WA^Mm 2 + \\A 1/2 v(t) II 2 < R 2 (36) 

with R and thus Cr independent of a and x. 

We consider a Lyapunov type function of the form 

*(*) = m+m, 

where 

E{t) = ^ (iK(t)H 2 + n\\A 1/2 w(t)\\ 2 ^ + Q(t) 

and 

<h(f) = rj(w, w t ) + fi(Kw, w), 

where 77 is a positive constant which will be chosen later and p = x+ija. By uniform dissipativity 
of the system (R,Sf’ x ) we can assume that that (36) holds with the same R as in (8) for all 
t > t *. 

One can see that there exists 0 < 7/0 < 1 and > 0 independent of (a; x) such that 

(3o[E 0 (t) + x\\K 1/2 w(t)\\ 2 - c R \\w(t)\\ 2 ] < 'L < foiEoft) + c R \\w(t)\\ 2 ] + / u||JL 1/2 7t;(t)|| 2 , 


for all t > f* and 0 < 77 < 770 , where 

E o(t) = \ (lK(t)|| 2 + n\\A 1 / 2 w (t) || 2 ) . 

Now on strong solutions we calculate the derivative 
dA> 

— = - ((D + aK)w t , w t ) - R(t) 

+ 77 [ 11 777 t 11 2 — (.Dwt , w) — v{Aw, w) — 2 x(Kw, w) — ( B{u) — B(v),w)] 
4 One can see that > x - inf spec(A'). Thus if K is not degenerate, then Sh —> +00 as k — > + 00 . 
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Since D is bounded from H 1 / 2 into H, we obtain that 

\(Dw t ,w)\ < e\\A l t 2 w\\ 2 + Ce~ l \\w t \\ 2 , Me > 0. 

Thus there exist 6,; > 0 independent of (a, x) such that 

^ < - [((£> + aK)w t ,w t ) - MIKII 2 ] + C R (\\u t \\ + \\v t ||)||A 1/2 u;|| 2 
- b 2 i][E 0 (t ) + x||A' 1//2 u;|| 2 ] + r]c R \\w(t)\\ 2 , 

Fixing a and taking s x large enough we obtain that 

^ + - C R (\\u t \\ 2 + \\v t \\ 2 )\\A 1/2 w\\ 2 <0, t> U, 

at 

for some u),C > 0. Using finiteness of the dissipativity integrals: 

(Kf + \\v t \\ 2 )dt < oo, 

we obtain (35). The equality 2l" ,x = 21 for the attractors follows from (35). □ 

If B is critical, but not satisfies the structural hypothesis in Assumption 3.7 we can still guarantee 
asymptotic exponential synchronization. However in this case we need additional condition that 
the damping parameter a is large enough and K is not degenerate. If B are globally Lipschitz 
we can even avoid the requirement of dissipativity of the system. 



Remark 3.15 The results similar to Theorems 3.12 and 3.14 can be also established for N 
coupled second order in time equations of the form 


u\ t + viAu 1 + D\uj + aK(u\ — u 2 ) + xK{v} — u 2 ) + B\(ii l ) = 0, 

u{ t + VjAui + Di‘u{ — aK(u{ +l — 2 u{ + u 2_1 ) — xK{vP +l — 2 u 2 + v?~ l ) + Bj(uP) = 0, 

j = 2,... , N — 1, 

u tt. + vnAu n + Dnu™ + aK(u^ — u ( V_1 ) + xK(u N — w w_1 ) + B^(u N ) = 0, 


This system can be reduced to (6) with 



(Di 

0 

0 

... 0 \ 


( l - 1 

0 

... o\ 



0 

D 2 

0 

... 0 


-1 2 

-1 

... 0 


V 0 = 

0 

0 

D 3 

... 0 

, K = 

0 -1 

2 

... 0 

K. 


\0 

0 

0 

... D N J 


\ o o 

0 

... 1/ 



Thus the general results of this section can be applied with the same hypotheses concerning the 
operators A, Di , K and H;. We note that the energy for for this N coupled model has the form 


N 

3=1 L 


(n«tii 2 + i/ jP 1 / 2 « j ii 2 ) + Uj(u j ) w r1/2 ( 


N-l 


(u j+1 — Vp) 


3 =1 
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In the ODE case (z^ = 0, K = id) synchronization for this model was considered in [24] with 
assumption that both a and x become large or even tend to infinity. Our approach allows us 
observe asymptotic synchronization for fixed a and in the limit x —> +oo (for identical subsystems 
it is sufficient to assume that x is large enough). The limiting (synchronized) regime is described 
by problem (27) with 


N N N 

j = 1 j= l i =1 

In the case of a plate with the Berger nonlinearity the same result was obtained in [32] with 
Dqj = dj ■ id, a = 0 and K = id. 

3.4 On synchronization by means of finite-dimensional coupling 

One can see from the argument given in Theorem 3.14 that asymptotic synchronization can be 
achieved even with finite-dimensional coupling operator. Indeed, the only condition we need is 

v{Aw, w) + x(Kw , w) > c||xo|| 2 , Vie £ H 1 / 2 , 

with appropriate c depending on the size of an absorbing ball. As it was already mentioned if K 
is a strictly positive operator, then the requirement s K > s* holds for large intensity parameter 
x. However it is not necessary to assume non-degeneracy of the operator K to guarantee large 
s x . For instance, if K = Hy is the orthoprojector onto Spanjefe : k = 1,2,..., N}, then 

N oo 

v(Aw,w)+x(Kw,w) > y>A fc + x)\(w,e k )\ 2 + v ^ A fc |(u;,e fc )| 2 

k —1 k=N-\-l 

N oo 

>(zzAi + \(w,e k )\ 2 + v\ N+ i ^ \( w , e k)\ 2 > min{z/Ai + x,v\ N+l }\\w\\ 2 . 

k =1 k=N-\-l 

Thus if x > v{ Atv+i — Ai), then we can guarantee large s H by an appropriate choice of N. 

This fact admits some generalization which based on an assumption that K is a “good” 
approximation (in some sense) for a strictly positive operator. 

Let V C H be separable Hilbert spaces and K, L be linear operators from V into H. Assume 
that L is a strictly positive on V, i.e., there exists > 0 such that 

( Lu,u ) > az,|M|fl-, u^V. 


We introduce the value 

e(L,K) = ey(L,K) = sup{\\Lu - Ku\\ h ■ \\u\\ v < 1}. 

In the case when L = id is the identity operator this value is known (see [2]) as the global 
approximation error in H arising in the approximation of elements v £ V by elements Kv. 

Now we take V = H 1//2 . It follows from the definition that 

||Lu - Ku\\ h < e(L, K)\\A l/2 u\\, u £ H 1/2 . 
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In this case we obtain 


u(Aw, w) + x(Kw , w) =z'||A 1 / 2 u>|| 2 + x(Lw, w ) + x((K — L)w, w) 

>zz||A 1//2 u;j| 2 + olx||w;|| 2 — xe(L, K)\\w\\\\A 1 ^ 2 w\ 
x 2 e 2 (L, K ) 


> a L x- 


Av 


w 


Thus according Theorem 3.14 under the condition 

x 2 e 2 (L, K) 


aL x 


4u 


> s* 


we have asymptotic exponential synchronization. The latter inequality is valid, when e 2 (L,I\) < 
2uaL>c~ 1 and x > 2s*af l for instance. So e(L,K) should be small and k large. We note that in 
the case L = id and K = Pn we have e(L, K ) = \ N +i and aL = 1- So the inequalities above can 
be realized for some choice x and N. 

Now we describe another situation where synchronization can be achieved with finite-dimen¬ 
sional coupling. For this we use interpolation operators related with a finite family C of linear 
continuous functionals {lj : j = 1,...., N} on H 1//2 . Following [9] (see also [19, 20]) we introduce 
the notion of completeness defect of a set C of linear functionals on H 1 / 2 with respect to H by 
the formula 


e c = ec{H 1/2 , H) = sup { || w \\ H - w € H 1 / 2 , l(w) = 0, Z G C, || w ||i/ 2 < l} • (37) 

It is clear that ec 1 > ec 2 provided Span£i C Span/l 2 and ec = 0 if and only if the class of 
functionals L is complete in H 1 ^ 2 , i.e, the property l{w) = 0 for all l G C implies w = 0. For 
further properties of the completeness defect we refer to [9, Chapter 5]. 

We define the class IZc of so-called interpolation operators which are related with the set of 
functionals given. We say that a operator K belongs to 7 Zc if it has the form 

N 

Kv = h (v)ifj, V v € H 1/2 , (38) 

3 =1 

where {ifj} is an arbitrary finite set of elements from H 1 / 2 . An operator K £ TZc is called 
Lagrange interpolation operator, if it has form (38) with {V’j} such that Ikif’j) = $kj- In the case 
of Lagrange operators we have that lj(u — Ku) = 0 and thus (37) yields 

\\u - Ku\\ < ec\\v\\ 1/2 , v£H 1/2 . 

Hence the smallness of the completeness defect is important requirement from point of view of 
synchronization. We refer to [9, Chapter 5] for properties of this characteristic and for description 
of sets of functionals with small ec- The simplest example is modes. In this case C = {lj(u) = 
(u,ej) : j = 1,2,... ,N}, where {e^} are eigenfunctions of the operator A. The operator Pc 
given by 

N 

P c y = ^2(e j ,v)e j , Vvi E i7 i/2 , 

3 = 1 
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_ 1 /2 

is the Lagrange interpolation operator. Moreover, e B = e{id , Pc) = • Thus the completeness 

defect (and the global approximation error) can be made small after an appropriate choice of N. 
This shows that the situation with K = Pjy can be included in the general framework. 

Unfortunately in the general case an interpolation operator of the form (38) is not symmetric 
and positive. Therefore formally we cannot apply the result on uniform dissipativity with K 
of the form (38). The situation requires a separate consideration and possibly another set of 
hypotheses concerning the model. We plan to provide with the corresponding analysis in future. 
Here we give only one particular result in this direction. 

We consider the following version of the equations (26) 

utt + vAu + Dut + xK(u — v) + B{u) = 0, (39a) 

vtt + vAv + Dvt + xK(y — u) + B(v) = 0, (39b) 

u(0) = n 0 , n t (0) = ui, u(0) = n 0 , n 4 (0) = v\. (39c) 

Theorem 3.16 Assume that A satisfies Assumption 2.1(i), B(u) is globally Lipschitz and sub- 
critical, i.e., 

3 CT 0 < 1/2 : \\B(ui) - B(u 2 )\ < L B \\A a °(ui - u 2 )\\, Vu* € H 1/2 . 

Let D : H 1 / 2 i->- H be strictly positive operator and K = Kc be a Lagrange interpolation operator 
for some family C of linear continuous functionals {lj : j = 1,, N} on H 1 ^ 2 . Then for every 
initial data Uq = (uo;uo) *= H l A and U\ = (ui;ui) £ H problem (39) has unique generalized 
solution U(t ) = ( u(t)]v(t )) and there exist x* > 0 and e*(x) such that for every x > x* and 
Lc < e*(x) the solution U(t) = (n(t);n(t)) is asymptotically synchronized, i.e., relation (35) holds 
with some ui > 0. 

Proof. This is globally Lipschitz case and therefore the well-posedness easily follows from [19, 
Theorem 1.5]. 

Following the same idea as in Theorem 3.14 we find that the difference w = u — v satisfies the 
equation 

w tt + uAw + Dw t + 2xw + G x xi u : v ) = 0> w(0) = w 0 , w t ( 0) = w\. 
where wq = uq — no, w\ = u\ — v\ and 

v) = —2 x{id — Kc){u — v) + B{u ) — B(v). 

We have the obvious estimate 

\\G K ,c(u,v)\\ < 2xec\\A 1 / 2 w\\ + Lb||„4 (T0 u;|| < 2{xe c + e>)||H 1 / 2 u;|| + C^IMI 
for all 5 > 0. This implies 

\{G h ,c{u,v),w)\ < (xe c + <5) 2 ||H 1/2 u;|| 2 + C' 5 ||n;|| 2 , V<5 > 0, 

and 

\{G K ,c(u,v),w t )\ < n\\w t \\ 2 + 2 ^ ec + ^ ||H 1/2 n;|| 2 + Cs,^\\w\\ 2 , Vp,6 > 0. 

/i 

Now as in the proof of Theorem 3.14 we can use Lyapunov type functional 
^(w,W t ) = 1 (jK|| 2 + u\\A 1/2 w\\ 2 ^ + x||u>i || 2 +7j(w,w t ), 
with x large, xec small and with an appropriate choice of p and 5. □ 
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4 Applications 

In this section we shortly outline possible application of the general results presented above. 

4.1 Plate models 

We first consider the plate models with coupling via elastic (Hooke type) links. Namely, we 
consider the following PDEs 

uu + 71 ut + A 2 it + xK(u — v) + <pi (u) = f\ in HcK 2 , (40a) 

v tt + 72 + A 2 u + xK{v - u) + tp 2 (u) = fi in SI C l 2 , (40b) 

with the hinged boundary conditions 

u = A u = 0, v = Av = 0 on <9P, (40c) 

where K is an operator which will be specified later. The nonlinear force term ip t (u) can take 
one of the following forms: 

• KirchhotT model: ip € Lipi oc (M) fulfills the condition 

liminf^ > -A?, 

|s|—>-oo S 

where Ai is the first eigenvalue of the Laplacian with the Dirichlet boundary conditions. 

• Von Karman model: <p{u) = [u,v(u) + To], where To is a given function in Tf 4 (P) and the 
bracket [u, v\ is given by 

[u, u] ' &X2 ^ ^X2 ^ ^XI^ ^ ’ ^X\X2^ ’ ^X\X2^ ’ 

The Airy stress function v(u) solves the following elliptic problem 

A 2 v(u) + [it, u] = 0 in P, — J' — = v(u) = 0 on <9P. 

on 

Von Karman equations are well known in nonlinear elasticity and describes nonlinear oscil¬ 
lations of a plate accounting for large displacements, see [20, 29] and the references therein. 

• Berger Model: In this case the feedback force has the form 

ip(u) = — k \Vu\ 2 dx — T A u, 

. Jn 

where k > 0 and T € R are parameters, for some details and references see [9, Chap.4]. 

In all these cases we have that H = T 2 (P) and 

Au = A 2 u, u € T>(A) = {«£ iT 4 (P) : u = A u = 0 on <9P} . 
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It is clear that A satisfies Assumption 2.1 (i) and T>(A 1//2 ) = H 2 (Q) n Hq(Q). 

The nonlinear force <p in the Kirchhoff model is subcritical with respect to the energy space 
(it is locally Lipschitz from H 1+S (Q) into Z, 2 (fi) for every 8 > 0). In contrast, the von Karman 
and Berger nonlinearities are critical (they are locally Lipschitz mappings from H 2 (Q) into L 2 (H) 
which are not compact on H 2 (Q)). Other requirements concerning the corresponding forcing 
terms Bi can be verified in the standard way. For details we refer to [14] for the Kirchhoff forces, 
to [19, Chapter 6] for the von Karman model and to [19, Chapter 7] for the case of Berger plates. 

The interaction operator K can be of the following forms K = id and K = —A, or even 
K = A a with 0 < cr < 1/2. In the purely Kirchhoff case with globally Lipschitz functions <pi we 
can also use the Lagrange interpolation operator with respect to nodes , i.e., with respect to the 
family of functionals 

lj(w) = w(xj), where Xj £ 12, j = 1,. .., N, 

with appropriate 5 choice of nodes Xj. This means that two Kirchhoff plates can be synchronized 
by finite number of point links. 

We note that in the case when both p>\ and ip 2 are Berger nonlinearities (possibly with different 
parameters) the results on synchronization with K = id can be found in [31], see also [32, 33]. 

We also mention some other plate models for which the abstract results established can be 
applied: 

• First of all we can consider the plates with other (self-adjoint) boundary conditions such 
as clamped and free and also combinations of them (for a discussion of these boundary 
conditions in the case of nonlinear plate models we refer to [20]). 

• The plate models with rotational inertia can be included in the framework presented. In¬ 
stead of (40) coupled dynamics in these models can be described by equations of the form 

(1 — yA )utt + Ml ~ yA)iq + A 2 u + xK(u — v) + (pi(u) = f\ in 12 C M 2 , (41a) 

(1 - yA )v tt + Ml — yA)?y + A 2 v + xK(v - u) + </? 2 (u) = f -2 in 11 C M 2 . (41b) 

Here y is positive. It is convenient to rewrite (41) as equations in H = Hq(Q.) (equipped 
with the inner product ( u,v)h = ((1 — yA)u, v) L 2 m)) in the form (26) with the operator A 
generated by the form a(u,v) = (Art, Au)L 2 (n) on LTq(H)- 

• We can also include into consideration the plates with Kirchhoff-Boussinesq forces of the 
form 

<p(u) = —div{ |Vu| 2 Vu} + a\u\ q u 

with a, q > 0, see [19, Chapter 7] and also [17, 21] concerning models with this force. 

4.2 Coupled wave equations 

In the case of coupled wave equations (14) the standard (critical) hypotheses concerning the 
source term ip £ C 2 (M) in the 3D case have the form: 

lirninf {(^j(s)s' 1 } > —Ai, \ip"(s)\ < C{ 1 + |-s|), sGi, 

|s|—HX> 

5 For details concerning smallness of the corresponding completeness defect we refer to [9, Chapter 5]. 
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where Ai is the first eigenvalue of the minus Laplacian with the Dirichlet boundary conditions, 
see, e.g., [19, Chapter 5] for details. 

We can also consider several versions of damped sine-Gordon equations. These are used 
to model the dynamics of Josephson junctions driven by a source of current (see, e.g., [39] for 
comments and references). For instance, we can consider the system 6 

utt + 7 ut ~ Art + /3u + x(u — v) + Asinu = f(x), (42a) 

vtt + 7 vt — Av + f3v + x(u — u) + Asinu = f(x), (42b) 


in a smooth domain D C and equipped with the Neumann boundary conditions 


du 

dn 


an 


= 0 , 


dv 

dn 


an 


= 0 . 


(42c) 


It is easy to see that in the case of the Dirichlet boundary conditions the general theory developed 
above can be applied. The same is true when /3 > 0. In the case f3 = 0 the situation is more 
complicated because the corresponding operator A is —A on the domain 

r c)i / 

V(A) = \ u£ H 2 (Q) : - = 0 on 
! on 


and thus becomes degenerate. So we concentrate on the case (3 = 0. 

It is convenient to introduce new variables 

u — v u + v 

w = - and z = -. 

2 2 

In these variables problem (42) with f3 = 0 can be written in the form 


(43) 


wtt + 7 w t 


Ztt + IZt - 


dw 

dn 



- Aw + 2xw + A sin w cos z = 0, 
Ax + A cos w sin z = f(x), 
dz | 


dn 


an 


= 0 . 


(44a) 

(44b) 

(44c) 


The main linear part in (44a) is not degenerate when n > 0. Therefore the same calculations as 
in the proof of Theorem 3.14 shows that there exists x* such that 


3t]> 0: \\w(t)\\ 2 Hl{n) + \\w t {t)\\ 2 <C B e rit , t> 0, 


when k > x* for all initial data from a bounded set B in I7 1 (D) x Lz^Q). This means that every 
trajectory is asymptotically synchronized. Moreover, it follows from the reduction principle (see 
[12, Section 2.3.3]) that the limiting (synchronized) dynamics is determined by the single equation 


ztt + ^zt — Ax + A sin x 


/(*)» 


dz 

dn 


an 


= 0 . 


The long-time dynamics of this equation is described in [39, Chapter 4], We also refer to [27, 28] 
for the studies of synchronization phenomena for (42) in the homogeneous (ODE) case. 

6 For simplicity we discuss a symmetric coupling of identical systems only. 


24 



Another coupled sine-Gordon systems of an interest is the following one 


u tt + 7 u t — Art + A sin(u - v) = fi{x), (45a) 

vtt + 7 v t - Av + A sin(u — u) = f- 2 (x), (45b) 

u\an = 0, u| 9 n = 0. (45c) 


Formally this model is out of the scope the theory developed above. However, using the ideas 
presented we can answer some questions concerning synchronized regimes. 

In variables w and z given by (43) we have equations 


Wtt + 7 wt — Aw + A sin 2 w = g(x), 

(46a) 

ztt + 7 zt - Az = h(x), 

(46b) 

w an = 0, z\qu = 0, 

(46c) 


where g(x) = (fi(x) — f2(x))/2 and h(x) = (fi(x) + f- 2 (x))/ 2 . One can see that 

|| 2 ,(t)l| 2 + ||V( 2 (t) - z,)f < C (|M0)|| 2 + IIVMO) - z,)|| 2 ) e-"‘, t > 0. 
for some C, ui > 0, where z* £ H l (Q) solves the Dirichlet problem 


-Az = h(x), z\ d n = 0. 


Problem (46a) equipped with the Dirichlet boundary conditions possesses a compact global at¬ 
tractor 21 in the space W = H ( \ (Q) x L^Q), see [39, Chapter 4], Hence 


(u{t) - zA 
u t (t) 
v(t) — z * 

V y t(t) J 


/ z{t) - z*\ 


/ w(t) \ 

ztit) 

+ 

W t (t) 

z{t) — z* 

-w(t) 

V z t(t) ) 


\-Wt(t)J 



^ £ 21 


as t-> +oo 


in the space W x W. Thus we observe some kind of shifted asymptotic anti-phase synchronization. 


A Appendix: Some facts from the theory of dynamical systems 

In this section we collect several definitions and general theorems from the theory of dissipative 
dynamical systems. For more complete presentation we refer to one of the monographs [3, 9, 23, 
39], 

Recall that a dynamical system (see, e.g., [9, 23, 39]) is a pair (X. St) of a complete metric 
space X and a family of continuous mappings St : X X, t > 0, satisfying the semigroup 
property: (i) St+ T = St o S T for any t, r > 0, and (ii) So is the identity operator. It is also 
assumed that 1 1-7 Stx is continuous mapping for every x £ X. 

A system (A, St) is said to be dissipative if it possesses a bounded absorbing set B. A closed 
set B C X is said to be absorbing for St if for any bounded set D C X there exists to(D) such 
that StD C B for all t > to{D). If the phase space A is a Banach space, then the radius of a ball 
containing an absorbing set is called a radius of dissipativity of the system. 

7 In the ODE case the synchronization phenomena in (45) was studied in [27, 28]. 
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A system (A, S t ) is said to be asymptotically smooth if for any closed bounded set B C X 
such that StB C B for all t > 0 there exists a compact set K. = IC(B) which uniformly attracts 
B, i.e., 

lim supjdist (Sty, 1C) : y G B} = 0. 

t— H-oo 

A system is called gradient if it possesses a strict Lyapunov function. The latter is defined as 
a continuous functional <b(x) on X satisfying (i) < <f>(x) for all t > 0 and x € X, and (ii) 

if <f>(x) = Q(Stx) for all t > 0, then x is a stationary point of St, i.e., St.x = x for all t > 0. 

A global attractor of a dynamical system (A, St) is a bounded closed set 21 C X which is 
invariant (i.e., Sfi t = 21) and uniformly attracts all other bounded sets: 

lim sup{dist(S'tj/, 21) : y G B } = 0 for any bounded set B in X. 

£—>•00 

It is known [3, 9, 23, 39] that the global attractor consists of all bounded full trajectories. We 
recall that a curve {y(t) : t G M} C X is said to be a full trajectory if S'ty(r) = y(t + r) for all 
rGl and t > 0. 

The standard criterion (see, e.g., [9, 23, 39]) for the existence of a global attractor is the 
following assertion. 

Theorem A.l Let (A, St) be a dissipative asymptotically smooth dynamical system on a complete 
metric space A. Then St possesses a unique compact global attractor 21 such that 

21 = u(B 0 ) = H U S - B ° 

£>0 t >£ 

for every bounded absorbing set Bq. 

For gradient systems it is also useful the following criterion of the existence of a global attractor 
(see, e.g., [35, Theorem 4.6]). 

Theorem A.2 Let (A, be an asymptotically smooth gradient system which has the property 
that for any bounded set B C A there exists r > 0 such that 7 T (B) = (J t>T St.B is bounded. If 
the set J\f of stationary points is bounded, then (A, S t ) has a compact global attractor 21 which 
coincides with the unstable set M + (AT) emanating from J\f, i.e., 21 = M + (AO- 

We recall (see, e.g., [3]) that the unstable set M + (A/") emanating from Af is a subset of A such 
that for each 2: G M + (A0 there exists a full trajectory {y(t) : t. G M} satisfying u( 0) = 2 and 
dist(y(f), J\f) 0 as t —>• —00. 

Remark A.3 We note that the hypothesis that 7 T (B) is bounded in Theorem A.2 can be changed 
in the following (additional) requirements concerning the Lyapunov function <h(x): (i) <h(x) is 
bounded from above on any bounded set, and (ii) the set = {1 G A : <3>(x) < R} is bounded 
for every R (see, e.g., [19, Corollary 2.29]). 

An important feature of a global attractor is its (fractal) dimension. We recall that the fractal 
dimension dinr^ M of a compact set M in a complete metric space A is defined (see, e.g., [39]) 
as 

dim A M = lim sup ————— , 

7 £ ^o hi(l/e) 
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where N(M,e ) is the minimal number of closed sets in X of diameter 2e needed to cover the 
set M. 

Now we state several facts related with asymptotically quasi-stable systems (for details we 
refer to the recent monograph [ 12 ] and the references therein). 

Let X and Y be reflexive Banach spaces; X is compactly embedded in Y. We endow the 
space H = X x Y with the norm 

\\ v\\ 2 h = IMIx + IMIy, y = («o;«i). 

Assume that ( H , St) is a dynamical system with the evolution operator of the form 

Sty = (u(t);u t (t)), y = (u Q ;ui) £ H, (47) 

where the function u(t) possesses the property 

u € c(R+,x) nc ,1 (M + ,y). 

A dynamical system (H, St) with an evolution operator of the form (47) is said to be asymptotically 
quasi-stable on a set B C H if there exist a compact seminorm yx{-) on the space X and 
nonnegative scalar functions aft), b(t), and c(t) on M + such that (i) aft) and eft) are locally 
bounded on [0, oo), (ii) bft) £ Li(R + ) possesses the property lim^oo b(t) = 0, and (iii) for every 
yi ,y 2 £ B and t > 0 the following relations 

\\S t yi - S t y 2 \\ 2 h < a(t) ■ \\yi - y 2 \\ 2 H 

and 

\\StVi - S t y 2 \\H < b(t) ■ \\yi - y 2 \\ 2 H + c(t) ■ sup [yxiu 1 ^) - u 2 (s))] 2 (48) 

0 <s<t 

hold. Here we denote Styi = {u l {t)\u\{t)), i = 1,2. The theory of quasi-stable systems was 
started with [15] and has been developed in [12, 19, 20, 22], The main outcome of this theory is 
the following assertion. 

Theorem A.4 (Global attractor) Assume that the system (H,St) is dissipative and asymp¬ 
totically quasi-stable on a bounded forward invariant absorbing set B in H. Then the system 
(. H,St ) possesses a compact global attractor 21 of finite fractal dimension. 

Another consequence of quasi-stability is the following assertion which states some regularity of 
the attractor and provides additional bounds for trajectories in it. 

Theorem A.5 (Regularity) Assume that the dynamical system ( H,St) possesses a compact 
global attractor 21 and is asymptotically quasi-stable on 21. Moreover, we assume that (48) holds 
with the function c(t) possessing the property c 0 0 = sup teR+ c(t) < oo. Then any full trajectory 
{(u(t); ut(t)) : t £ M} that belongs to the global attractor enjoys the following regidarity properties 

Ut £ LoofM.', X) C\ C(M.]Y), u tt £L 00 (R-,Y). 

Moreover, there exists R > 0 such that 

H^tWIIx + ll^ttWIly ^ R 2 1 t £ M, 

where R depends on the constant c QQ , on the seminorm fix, and also on the embedding properties 
of X into Y. 
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In the study of synchronization phenomena we deal with stability of attractors with respect 
to interaction intensity. Variations of the intensity parameter is treated as a perturbation of 
the dynamical system. To describe these variations at the abstract level we consider a family 
of dynamical systems (X,S^) with the same phase space X and with evolution operators S A 
depending on a parameter A from a complete metric space A. The following assertion is proved 
by Kapitansky and Kostin [26] (see also [3] and [23] for related results). 

Theorem A.6 (Upper semicontinuity) Assume that a dynamical system (X,Sf) in a com¬ 
plete metric space X possesses a compact global attractor 2l A for every A € A. Assume that the 
following conditions hold. 

(i) There exists a compact K C X such that 21 A C K. 

(ii) If Afc -* Ao, Xk —>• xq and Xk € 2l Afc , then 

S^ k Xk —> S^°xq for some r > 0. (49) 

Then the family {2l A } of attractors is upper semicontinuous at the point Ao/ that is, 

dx |2l A 121 A °| = sup jdistx(u2l A °) : x G 2l A | —> 0 as A -» Ao- 

Moreover, if (49) holds for every t > 0, then the upper limit 2l(Ao, A) of the attractors 2l A at Ao 
defined by the formida 

2l(A 0 ,A) = run* : A e A, 0 < dist(A, Ao) < 5} 

5>o 

is a nonempty compact strictly invariant set lying in the attractor 2t A ° and possessing the property 

dx "12l A | 2l(Ao, A) ^ —> 0 as A —> Ao• 
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